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O : Abstract 

We consider two three-dimensional isotropic harmonic oscillators interacting with the quantum 
electromagnetic field in the Coulomb gauge and within dipole approximation. Using a Bogoliubov- 
like transformation, we can obtain transformed operators such that the Hamiltonian of the system, 

cr 

when expressed in terms of these operators, assumes a diagonal form. We are also able to obtain 
^ . an expression for the energy shift of the ground state, which is valid at all orders in the coupling 

h : 

constant. Prom this energy shift the nonperturbative Casimir-Polder potential energy between the 
two oscillators can be obtained. When approximated to the fourth order in the electric charge, 
the well-known expression of the Casimir-Polder potential in terms of the polarizabilities of the 
oscillators is recovered. 
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I. INTRODUCTION 



The Casimir- Polder potential is a long range interaction between neutral atoms/molecules 
in vacuo This potential arises from the interaction of the two atoms with the common 
radiation field. It is also a result of the quantum field fluctuations in the vacuum state, and it 
is one of the few unambiguous manifestations of the quantum nature of the electromagnetic 
field Q. 

The Casimir-Polder potential is usually obtained by fourth-order perturbation theory Q]. 
Higher-order corrections exist, as well as non-additive components in the case of three or 
more atoms \m. More recently Casimir-Polder forces for excited atoms have been considered 
B, eve n if the precise meanlng of a dressed excited state is not clearR fl g as well 
as time-dependent Casimir-Polder energies for partially dressed atoms |9j or atoms in an 
excited-state [loj | . 

The Casimir-Polder potential can be expressed in terms of the dynamical polarizabilities 
of the two atoms/molecules, and this makes clear that a similar potential is expected when- 
ever we consider two neutral systems with a discrete spectrum interacting with a quantum 
field (not necessarily in the framework of quantum electrodynamics, but also in condensed 
matter physics, for example). 

In this paper we consider two harmonic oscillators interacting with the electromagnetic 
radiation field. This system has the advantage that exact solutions can be obtained by a 
Bogoliubov-like transformation, as we will show in this paper. It is also a good representation 
for atomic systems; approximating atoms by harmonic oscillators is a common procedure, 
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due to the resulting simplification of calculations (for a few examples, see 
An harmonic oscillator interacting with a one-dimensional scalar field has been recently 
considered in order to study phenomena relevant to quantum optics, and exactly solved by 
a Bogoliubov transformation 0, ^| ; similar methods have been used also for the description 
of unstable states in relativistic field theory An exact solution for two two- level atoms 
interacting with a one- dimensional scalar field, in the one excitation sector, has been recently 
obtained for the description of decaying states in the rotating wave approximation Q| . The 
rotating wave approximation, however, prevents application of these results to the Casimir- 
Polder potential, where the virtual processes play an essential role. 

We are mainly interested in the study of the Casimir-Polder potential. We therefore con- 
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sider two spatially separated isotropic three-dimensional harmonic oscillators, (thus having 
a spherical symmetry), interacting with the quantum electromagnetic field. Counterrotating 
terms and the related virtual processes are included in the Hamiltonian. By a Bogoliubov- 
like transformation we can obtain an exact diagonalization of the Hamiltonian describing 
this system. The Hamiltonian of the system, in terms of the transformed operators, has the 
form of a free field Hamiltonian. Moreover, using the inverse transformations we are also 
able to obtain the exact energy shift of the ground state of the system, due to the interaction 
between the oscillators and the radiation field. This result is exact within the model consid- 
ered and valid at all orders in the coupling constant. The part of the energy shift depending 
on the distance between the two oscillators yields the Casimir-Polder potential between the 
ground-state oscillators, both in the so-called near and far zones. Known results are repro- 
duced by approximation to the fourth-order in the electric charge. Due to the generality of 
our model, and the fact that the validity of our results is not limited by perturbation theory, 
the results obtained in this paper can be also used to obtain Casimir-Polder-like potentials 
in all cases in which the coupling constant is not small. We also expect that our nonper- 
turbative approach can be useful in dealing with Casimir-Polder forces for atoms in excited 
states, being able to overcome the conceptual difficulties related to the vanishing energy 
denominators occurring in the perturbative expansions when excited states are considered 

1 E21. 

In Section |H] we introduce the Hamiltonian describing our system and the transformation 
used to obtain the new operators; in terms of the new operators, the Hamiltonian assumes a 
diagonal form. The explicit relation between the new and old operators is also obtained, as 
well as the inverse relations. In section HTT1 we obtain the energy shift of the ground state of 
the interacting system, from which the Casimir-Polder potential between the two isotropic 
harmonic oscillators is obtained. The physical meaning of our results is discussed, as well 
as possible applications to Casimir-Polder forces for excited atoms. 



II. HAMILTONIAN AND BOGOLIUBOV TRANSFORMATION 

We consider two isotropic three-dimensional harmonic oscillators with frequency uq = cko 
interacting with the electromagnetic field in the Coulomb gauge. In the multipolar coupling 
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scheme the Hamiltonian of this system is 

#=E E k a\ e a te + Y,kBi 3 B kj +Hl t + Hi t (1) 

i=l,2£=x,y,z kj 

(we use units such that h = 1 and c = 1). The index i labels the two oscillators (i = 1,2), 
the index I the space direction of the oscillator (£ = x,y,z), al^aa are the creations and 
annihilation operators of the two three-dimensional oscillators and B^pB-^j for the field (j 
is the polarization index), k —\ k |. The creation and annihilation operators in satisfy 
the usual bosonic commutation rules 



a i£? a \'£> 



hv&UJ i ^kji-BjL, — <W%'- (2) 



H} nt and Hf nt are the interaction Hamiltonians of the oscillators 1 and 2 with the field, 
respectively. Within the dipole approximation, they are given by 

hL= E E/^ + aL) K-^y (3) 

£=x,y,z kj 

hL= E E/kK^ + 4) Ke fa -4ie- ,k j (4) 

£=x,y,z kj 

where the atom 1 is supposed to be located at the origin and the atom 2 at r, and the 
coupling constants are 



\x it is the transition dipole moment of the atom % in the £ direction and V is the quantization 
volume; the dipole moments are supposed real. 

In order to diagonalize the Hamiltonian we look for a linear transformation of the 
operators, similarly to the case of just one one- dimensional oscillator interacting with a 
scalar field Thus we write the new operators frkj^k? in terms of the old operators in 
the following form 

bh = E E (*&4 + rgotf) + E (^Xy + (6) 

i=l,2 £=x,y,z k'j' 

and its Hermitian conjugate for the annihilation operators b^j. The coefficients 
^kjJ r kj'^k/ ' -^k/ are chosen in such a way that the "new" operators satisfy the free-field 
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commutation relations 



[H,b kj 



kbl 



hi 
kb kj 



(7) 



Substitution of JBJ) into (J7|) yields the following set of coupled equations for the coefficients 



(k - fco)*g 


= J2fk'j' 




+0 




(8) 


(k + k Q )r% 


frpk'j 

V^j 






(9) 


(k - k Q )t% 


= J2fk/j' 
k^' 

= J2fk/j' 
k'j> 


frpk'j 

V^j 


e* k ' r + it*/e- ik '- r ) 




(10) 


(k + k Q )r% 


frpk'j 

V^j 


e ik '- r + it! k /e- ik ' r ) 




(11) 


(k - k')T«/ 


= J2fk'j' 




-«)+E/k 2 ^ k ' r (^ 


"«) 


(12) 


(k + k')R«f 


= E fk'j' 


(r U 






(13) 



This set of coupled equations can be solved exactly. The solution can be expressed in terms 
of the resolvent 

<?.k' .. .. 

(14) 



(CC*)) -1 = kl -k 2 - 2k T^^fk'ffk'j' 
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and the r-dependent function 



a e (k,r) = -2k J2 



gik'r „— jk'-r 



1 f U r2£ 
7 Jk'j'Jk'j' 



3 V 



(15) 



For isotropic oscillators, as we assume, the resolvent (|14j) does not depend on t. Explicit 
expressions of (fT4*|) and ()15|) in the continuous limit are given in the Appendix [X] In this 
limit, the two functions (|T4*j) and (fTo^) have poles for k = k', and must be extended to the 
complex plane 



G { ^\k) = Gi(k±ie) 
a^\k,r) = <jf(k±ie,r) 



(16) 
(17) 



However, whenever unnecessary, for simplicity of notations we shall omit the apix relative 
to the choice of the analytic continuation of these functions. For convenience of notations 
the two functions defined in (|T4*|) and (fTB^) can be expressed as the elements of a symmetric 



2x2 matrix, with (Gi) 1 and (G2) 1 as diagonal elements and at as off-diagonal elements, 
apart from a proportional factor. So we define the matrix 

GiG 2 (°2 l -at 
1 — afGiG 2 \ —a t Gi 

We also define the following two-dimensional column vectors 



* = U1U (i9) 

*-(*)' *-Cs) ^ ^ (20) 

The following relations exist among and tj^, and among i?^ and T k - J 

= 1^* < 21 » 

%f = Irl^ < 22) 

The solution for the coefficients defined in © can be expressed in the compact form 

t£j = -(* + ko)M,f<i (23) 

r£j = "(t-MM,^ (24) 

j£ ¥ = M* + J3p (4,)' M ( f^ (25) 

< / = -F^( f ^) ,M '^ < 26 > 
The correct commutation relation for the transformed operators 



SjfSw (27) 



can be proved by expressing b^j and in terms of the "old" operators by means of eq. (JHJ) 
and eqs. (|23I24I25I26|) . Equations ()A7|) and (|A8|) of Appendix |X] have been used in order to 
obtain the commutation relations ([27)1 . 

The inverse relations, expressing old operators in terms of new operators, can be also 
obtained in the form 

,t _ V" 



(A!? bl + ^ b{ 3 ) (28) 

Bl 3 = E (/#' + M (29) 
k'i' 
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and the Hermitian conjugate relations for an and B^j. The coefficients appearing in eqs. 
(128129)) can be obtained with the same procedure used for the "direct" relations. The result 
for the A coefficients is 

^ - r^lmkm - fi^r)G l(k) ) (3D 

The other coefficients may be expressed in terms of the A coefficients as 

2 = MO* P2) 



T " k + h, - ■ 

# = FTI W (34) 



III. THE GROUND-STATE ENERGY 

We now evaluate the energy shift of the ground-state energy due to the interaction of the 
two harmonic oscillators with the electromagnetic field. The ground-state energy has terms 
which depend on the distance r between the two oscillators; this yields a Casimir-Polder 
potential between the two oscillators. 

In order to evaluate the energy of the ground state we substitute the inverse relations 
()3(JI31l32l33l34j) in the Hamiltonian (|T|). After this substitution, the Hamiltonian in terms 
of the new operators assumes a diagonal form 

H = J2kbl J b kj + E (35) 

where Eq is the ground-state energy. The part of Eq depending from the interatomic distance 
r gives the Casimir-Polder potential energy between the oscillators. After extensive algebraic 
calculations, and use of eqs. (|A9|) and (|A10|) of Appendix |XJ we obtain 



(it should be note that the A functions in (jHUj) depend on r) . This expression is valid at any 
order in the electric charge. In the continuous limit, using the results of Appendix El we 
can obtain an explicit expression of the shift of the ground state energy E . By retaining 
only the r dependent parts of eq. (J36j) up to the 4 th order in the charge, approximating to 
the far zone (k, k! -C /cq) |2fl, the following expression for the interatomic potential Vcp{r) is 



obtained 



VWr) = ^(2/ 1 + /,)i (37) 



where I\ and I 2 stand for the following integrals 



/•OO POO 

h = dx dy 
Jo Jo 



1 A 

x y 
x 2 — y 2 



OO POO r^nA 



3 POO 

dx I dy- 
Jo 



x"y- 



[x + y) 



^)(m- m )+j^)( 3JM -Mv) 



(38) 



(39) 



Once the above integrals are computed by standard methods as l\ = — ^ and I 2 = ^r, 



the final result is 



VcW - (40) 



where ai (i = 1, 2) is the static polarizability of the i-th three-dimensional isotropic harmonic 
oscillator, that is 

a . = \ ^ = ^t (41) 

O „ kn fcn 

t=x,y,z u u 

The quantity (J4*U|) is the well-known expression of the Casimir-Polder potential in the far 
zone and for isotropic molecules . Therefore, in the appropriate limits, the exact result 
(jnnj) of the ground-state energy E yields the usual Casimir-Polder potential energy and this 
confirms the validity of our approach. 

IV. CONCLUSIONS 

We have considered a system of two isotropic three-dimensional harmonic oscillators in- 
teracting with the electromagnetic radiation field, described in the Coulomb gauge and in 
the dipole approximation. Through an appropriate transformation we are able to diagonalize 
the Hamiltonian of this system and obtain exact solutions. The diagonalized Hamiltonian 
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is an infinite set of harmonic oscillators with the same spectrum of the free electromagnetic 
field, plus an energy shift. We have obtained an expression of this energy shift valid at 
all orders in the coupling constant, which contains terms depending from the distance r 
between the two three-dimensional oscillators. These r-dependent terms yield a potential 
energy between the two oscillators in the ground state of the interacting system, i.e. the 
Casimir-Polder interaction. Our results allow to obtain this energy exactly, at any order 
in the coupling constant. Approximation to the fourth order in the coupling constant re- 
produces known results of the Casimir-Polder long-range interaction in terms of the static 
polarizabilities of the oscillators (in the far zone). We plan to discuss in a future publica- 
tion the extension of the results obtained in this paper to the case of the Casimir-Polder 
potential for excited states. In this case, we should be able to obtain the potential without 
the well-known difficulties related to vanishing energy denominators that appear when the 
stationary perturbation theory is applied to excited states. 

Finally, we wish to stress that our method is quite general and not limited to systems 
interacting with the electromagnetic field. Our results can be easily extended to other 
matter-field interacting systems, for example the interaction with phonons. Thus the method 
presented in this paper could be particularly useful for various physical systems with a strong 
coupling constant. 

APPENDIX A: THE FUNCTIONS G^k) AND a(k,r) IN THE CONTINUOUS 
LIMIT 

In this Appendix we derive explicit expressions for the functions defined in (|14j) and (|15|) . 
in the case of isotropic oscillators, as well as quantities useful for the explicit evaluation of 
the commutation relations (}2T)) and of energy shift E given in Sec. III. The hypotheses 
of isotropic oscillators allows significant simplifications of the final equations. In equation 
(fTlj) the summation over the polarization index f can be done using Y!,j'(^k'j')m(^k.'j')n — 
Smn — k' m k' n . Also, in the continuous limit 




(Al) 



Thus, after polarization sum and angular integration we obtain 




(A2) 
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where we have used that for symmetrical oscillators the strength of the dipole moment in 
the three spatial directions is the same, here indicated with /ij (i = 1,2). 

We can proceed similarly for the quantity appearing in (|15j) . After summation over 
polarizations and angular integration, we obtain 



(A3) 



The function he(k'r) is defined by 



he(k'r) 



jo(k'r) - ±-ji{k'r) for £ = x,y 



(A4) 



p^'i(fc'r) for^ = 2; 

where the axis z has been taken along the direction of r, and jo,ji are spherical Bessel's 
functions. 

Thus, in the case of isotropic oscillators the functions f)14j) and (fT3j) take the form 

4 . o r „, 2k' 4 



mk)))- 1 = k 2 Q - k 2 - -^kofij dk'y 2 



k 2 

2 r 2k' 4 
ae(k,r) = --ko/j,^ \ dk' r^h(k'r) 



(A5) 
(A6) 



7i J k — k 2 

Equation ()A5j) contains an ultraviolet divergence that, however, is inessential for our pur- 
poses. flA6|) is convergent due to the oscillating behaviour of hi(k'r). 

By using the definitions (fT4^) and (fT3j) of (G^A;)) -1 and ae(k,r), the following identities, 
which have been used to obtain the commutation relations ()27|) and to evaluate the energy 
shift E , hold 



kj 



(Gi(k')y 1 - (Giik) 



1 

2k /,■' - k 

cr £ (k',r) - ae(k,r) 
2k (k - k') 



-i 



+ (k' + k) 



(A7) 
(A8) 



ki 

ki 



A' 
~2 



k >z + k s + k' (G l (k')y 1 + k(G l Ck) 
2k (k' + k) 

k'ai(k' , r) + kae(k, r) 



2k (k> + k) 



(A9) 
(A10) 
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where F^(k) and F^(k) are given by 



F (D (k) = 1 1 = (_}: L_ 1 fA ll) 

k '' kK ' (k'-k)(k-k) (k + k)(k' + k) k-k'\k 2 -k 2 k 2 -k' 2 ) v ; 

F m (k) = 1 , 1 = 2 _ ( k ' , _J_\ (A12) 

k '> kK ' (k'-k)(k + k) (k-k)(k' + k) ~k + k'\k 2 -k' 2 k 2 -~k 2 J K ' 
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